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Adsorbed noble gas atoms donate (on the average) a frac- 
tion of an electronic charge to the substrate metal. The effect 
has been experimentally observed as an adsorptive change in 
the electronic work function. The connection between the ef- 
fective net atomic charge and the binding energy of the atom 
to the metal is theoretically explored. 

PACS: 73.30.+y, 34.50.Dy, 68.10.Gw 



I. INTRODUCTION 

The problem of measuring and understanding the 
binding energy of noble gas atoms on metal substrates 
has always been of considerable interest The 
atomic binding is that energy released when an atom 
from the vapor sticks to the surface. Several studies 
have related atomic binding to surface charge distribu- 
tions [l|,|7j-p^|| . Induced surface electronic dipole moments 
near the adsorbed atom have been of particular interest. 
In the work which follows, the intimate relationship be- 
tween binding energy and induced atomic charge will be 
considered in theoretical detail. 

As the atom is lowered onto the surface and becomes 
adsorbed, the atomic dipole tends to be oriented with the 
positive side of the dipole pointing away from the metal. 
In reality, the negative side of the atomic dipole is better 
thought to be on average an electron (negative) charge 
—Z e ff\e\ donated to the metal. This leaves the atom with 
a positive net mean charge +Z e ff\e\. The physical situa- 
tion is pictured in FIG.l below. Even in a situation often 
regarded as physisorption, one does not expect a noble 
gas atom to remain in perfect charge neutrality. When 
the atom is adsorbed on the substrate, the negative end 
of the atomic dipole moment neutralizes the positive end 
of the image dipole moment. This leaves a mean net pos- 
itive charge on the atom and a mean negative electronic 
charge deposited in the metal. From an experimental 
viewpoint, the positive nature of the mean charge on an 
adsorbed noble gas atom is equivalent to the donation 
of a negative charge to the metal. The physical effect is 



made manifest by the diminution of the 

electronic work function as the hrst monolayer of atoms 
is deposited on the metallic substrate surface. This effect 
is quite large and is observed for all combinations of gases 
and metals . The magnitude of the reduction of the 
work function is proportional to the amount of adsorbate 
and for coverage up to a monolayer Jlf| . 
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FIG. 1. Shown is the surface between the vacuum and the 
metal. When an atom is far above the surface between the 
vacuum and the metal the dipole fluctuations are correlated 
with the image dipole, as shown to the left. As the atom is 
lowered onto the surface, the atomic dipole is pointed upwards 
as shown on the right. An electron is then donated to the 
metal with probability giving the effective charge strength 

In Sees. II and III the long ranged part of Van der 
Waals interaction between a noble gas atom and a metal 
will be reviewed. The height of the atom, over and above 
the metal surface, is considered to be large compared 
with the atomic size. The interaction is between a quan- 
tum fluctuating dipole and its correlated image. Both 
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the binding energy and the probability of electronic ex- 
citation will be computed numerically in Sec. IV. In Sec. 
V, non-perturbative expressions for the binding energy 
and the probability of electronic excitation are derived. 
In Sec. VI it is shown that the stronger the binding en- 
ergy the larger the net atomic adsorptive charges. The 
physical importance of the adsorptive charge on a noble 
gas atom is discussed in the concluding Sec. VII. 



II. PERTURBATION THEORY 

Suppose a Hamiltonian of the form 

H = H + V, (1) 

in which 

Hiji n = E n ^ n , (2) 

represents an unperturbed energy eigenvalue problem, 
and 



(3) 



represents the perturbed energy eigenvalue problem. 
Without loss of generality, one may assume (for the un- 
perturbed ground state wave function) that 



(Vo,V^) = 0. 



(4) 



Two quantities of importance are the exact ground 
state energy £q, and the probability that the interaction 
V introduces an excitation in the unperturbed quantum 
states 



p = ^i(^„,* )i 2 = i-i(^o,*o)r 

n#0 



(5) 



In lowest order perturbation theory, one then finds the 
energy shift H 



U=(£ Q -E Q ) = -J2 



(E n — E ) 



(0) 



The ground state wave function in first order perturba- 
tion theory 20 1 is given by 



{E n — Eci) 



(J) 



The first order wave equation can be used to compute 
the excitation probability P to second order in the per- 
turbation potential; it is 



P 



E 

n#0 



On, Vlp ) 



{E n — Eq) 



(8) 



Let us apply these ideas to the Casimir effect, i.e. the 
attractive force on an atom located at height h above 
a metallic substrate. The final energy is well known. 
However, the excitation probability of the atom will also 
be calculated. 



III. THE VAN DER WAALS FORCE 

Let the unperturbed problem consist of an isolated 
atom and an isolated metal 



H = H n 



Hr, 



;etal • 



(9) 



With n = {j, a) as a double index referring to the atomic 
state j and the metallic substrate state a, we may write 



i>n = vf om vC etaZ , 



with an unperturbed energy 
E 



j-i a torn , rpmetal 

&j + K 



(10) 
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Further, let the interaction between the atom and 
metal be of the dipole form 



V = -/j,-~E 



(12) 



where /i represents the atomic electric dipole moment 
operator, and the electric field E is produced by the 
metal. The second order perturbation interaction en- 
ergy between the atom and the metal is usually called 
the Casimir effect and is given by 



U c 



E 



(aj)#(0,0) 



(hujjo + hu' a0 ) ' 



where 



and 



» j0 = &f° m ,^ tom ), 



E/.i,metal Tij,metah 
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The atomic Bohr frequencies are defined as 
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(Ef° 
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and the metallic Bohr frequencies are defined as 



E 



metal 



E 



metal \ 



(13) 

(14) 
(15) 
(16) 
(17) 



For an isotropic atom, one may define a ground state 
dipole moment quantum noise spectral function 

5 M = I E I «° m ' ^o t0m ) f 6 w io) , (18) 

Similarly, one may define the electric field quantum zero 
point fluctuations due to the metal 

- E I hPa etal ^Z tal )\ 2 ^ (U/ - W' a0 ) . (19) 
a#0 



One may now compute the Casimir energy j21|-j23| Uc in 
terms of these spectral functions; i.e. 



2 



Finally, one may find the spectral functions by employing 
the zero point quantum fluctuation response theorems 
II- 

For example, if the atomic polarization response to 
an external electric field at complex frequency £ (with 
3m(C) > 0), 



S (jj.) = a(C)5E ea;t , 



(21) 



defines the ground state atomic polarizability ct(£) for 
a spherical atom, then the fluctuation response theorem 
theorem asserts 



(22) 



S(u) = - 3?m a{uj + i0+). 



The fluctuation spectral function for the electrostatic 
field E = — V0 produced by the metal is a bit more 
subtle. 

If an external charge density at complex frequency £ 
produces an electrostatic potential according to the rule 



* (#r)) = / <?(r,r',C)<W(r')dV, 



(23) 



then the spectral function for zero point electrostatic po- 
tential fluctuations 



5 (r, r',ui) 



dt 



cos(^)5Re (0 |A<£(r, t)A(f>(r', 0)| 0) ( — J (24) 



obeys the fluctuation response theorem in the form 

S^r, r', u) = - (^j 3me(r, r', uj + »0+). (25) 

For the problem at hand, suppose that the metal is lo- 
cated in the half-space z < 0; e.g. the substrate surface is 
the z — x-y plane. If both r and r' are in the vacuum 
(i.e. z > and z' > 0), then the "method of images" 
yields the Greens function 



S(r,r',C) 



1 



(26) 



Jr — r'| ) \ |r — r,- 
where 

r' = (x',y',z% (27) 
has a corresponding "image" position 

v[ = {x\y',-z'), (28) 

and 

e(0 1 



v(0 



<0 + i) 



(29) 



The dielectric response function for the metal e(Q deter- 
mines the conductivity cr(£) via 



KC) = i + (^) 



(30) 



Thus 



^(r, r', u) = ( ^3^1 ) 3m + i0+ )- ( 31 ) 



With 



R= (0,0, h) 



(32) 



denoting the position of the atom at a height h above the 
substrate surface, the electric field zero point fluctuations 
of the electric field at the atom may be computed via 

Se(w) = lim lim 

r^R r'^R 

/ d 2 d 2 d 2 \ 

( v 9^ + W + ^J^ (r ' r>) - (33) 



The differentiation is tedious but direct. It yields 



5e(w) = (27 Tfi 



(34) 



Substituting Eqs.(22) and (34) into Eq.(20) yields the 
zero-point fluctuation Casimir effect potential 



U C = - 



C a 



where 



Catom — 



h 3 



2tt 2 



(35) 



JQ 



3m a[ui + i0 + )Sm rj(u)' + i0 + ) 

LO + LO' 



dtvduj' . 



(36) 



The Casimir form of the Van der Waals potential between 
an atom and a conducting surface is well known. The 
purpose for reviewing the Van der Waals result is that 
now we may also calculate the probability that the atom 
at a height h is in an excited state. 

The excited state probability in Eq.(5), when evaluated 
to the lowest order perturbation theory in Eq.(8), yields 
the Casimir excited state probability Pq ■ In the Casimir 
perturbation theory, the probability for the atom to be in 
an excited state contains an extra energy denominator. 
The final result is in the simple form 



Pc 



/Va_ 



h 3 



(37) 
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The parameter v at0 m is a rough measure of the effective 
atomic volume; It is 



1 



Vatom — 1 ' X 



3m a{ui + ?0 + )3m r)(u>' + i0 + ) 



dtvduj' . 
(38) 



The excited state probability in Eq.(37) is a new result. 
The noble gas atom above the metal, is not electronically 
inert. The atom in the vacuum is in the ground state, As 
the atom is lowered toward the metal surface, the atomic 
electronic state becomes excited with probability P > 0. 
For heights h such that h 3 >> v a tom we have P s» Pc as 
in Eq.(37). Numerical results for the probability P c then 
follow (below) in a similar manner to numerical results 
for Uc which have been previously computed by other 
workers. 



IV. NUMERICAL EVALUATIONS 

For the purpose of numerical evaluations of parameters 
associated with the Van der Waals interaction, we follow 
Rauber, Klein, Cole and Bruch |f25[ ; They choose a single 
pole approximation for both a(£) and 77(C)- The pole in 
a(C) is parameterized by an atomic frequency u> a , while 
the pole in 77(C) is parameterized by a plasma frequency 

n,. 

In detail, the approximations read 



, RKCB 



(0 



along with 



RKCB 



n 



(0 



u„ — L 



Eqs.(36), (38), (39) and (40) yield 



RKCB _ ft J CtQgoUgSls 



and 



C. 



RKCB 
atom 



8\" 



fi. 



1 f a Q goUj a n s 



(39) 



(40) 



(41) 



(42) 



The values of the atomicparameters ao and uj a have 
been previously tabulated f2q| . and are listed in Table I. 
The values of go and f2 s , (for several metals) have also 
been previously tabulated j25) and are here listed in Table 
II. In Table III, we have computed the values of C a tom 
and v atom - 

The probability of atomic excitation in the Casimir 
regime Pc = {v a tom/h 3 ) can be calculated in the limit of 



large h. When the height of the atom obeys h 3 » v atom , 
the perturbation theory is reliable. However, for atoms 
adsorbed on a submonolayer of film the height h is not 
large. Thus a non-perturbative method is required. We 
now turn our attention to this more detailed treatment. 

TABLE I. Values of qq and hw a for Noble gas atoms. 



Atom 


(qo/A 3 ) 




(hcoa/eV) 


He 


0.205 




27.645 


Ne 


0.396 




32.734 


Ar 


1.642 




18.971 


Kr 


2.487 




16.478 


Ac 


4.018 




1 A O A 

14. 61 


TABLE II. 


Values of go and hn„ for some 


metals. 


Atom 


(an) 




(hQ 3 /eV) 


Cu 


0.857 




17.74 


Ag 


0.812 




21.768 


Au 


0.840 




24.162 


Al 


0.976 




12.87 


Pd 


0.785 




17.115 


Gr 


0.619 




18.149 


TABLE III. 


Catom &nd Vatom 


for noble gas 


atoms. 


Gas-Metal 


{Catom/ -A^6 


V) 


(Vatom/A 3 ) 


Xe-Cu 


3.413 




0.106 


Xe-Ag 


3.525 




0.098 


Xe-Au 


3.797 




0.099 


Kr-Cu 


2.276 




0.067 


Kr-Ag 


2.367 




0.062 


Kr-Au 


2.558 




0.063 


Ar-Cu 


1.612 




0.044 


Ar-Ag 


1.690 




0.041 


Ar-Au 


1.832 




0.043 


Ne-Cu 


0.488 




0.096 


Ne-Ag 


0.525 




0.096 


Ne-Au 


0.578 




0.010 


He-Cu 


0.237 




0.005 


He-Ag 


0.253 




0.005 


He-Au 


0.277 




0.005 
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V. RIGOROUS RESULTS 

The above considerations are true for atoms above a 
metal in the perturbative limit h — > oo. When the atom 
is adsorbed on the metal surface, a non-perturbative 
viewpoint must be invoked ^f^7|. For example, con- 
sider the ground state matrix element analytic in the 
upper half complex energy plane 3m z > 0, 



G(z) = 



1 



z-H. 

where -00 is the unperturbed ground state, 

Hi/jq = E ip , 
and the binding energy of the adsorbed atom, 

U = -B, B>0, 
is determined by the full electronic ground state 
H^fo = (H+ V)V = {E + U)V . 
One may write Eq.(43) in the exact form 

1 



G{z) = 



E - E(z) 



(43) 

(44) 
(45) 
(46) 

(47) 



where the self energy part determines the atomic 

binding energy 



U = T,(E + U). 



(48) 



G{z) has a simple pole at the exact ground state energy 
z = E + U = E — B. The residue at the pole is the 
transition probability ipo — > ^oJ i-e. 



l(*o,^)| s 



1 



1 - E>(Eo + U) 



(49) 



where £'(2) = dE(z)/dz. The exact expression for the 
self energy part reads 



i 



(50) 



where TC = PTLP and where P projects into the subspace 
normal to the unperturbed ground state. For example, 
from Eq.(49) and the definition of P one finds the mean 
value P — (^o, P^o) given by 



P = l-|0&o,^o)| 2 = 



1 - E'(£7 + U) 



(51) 



The details of the mathematical derivation of Eqs.(47)- 
(51) are given in the appendix. 

Suppose that the atom (located on the metal surface) 
starts out in the unperturbed state tpQ. The atom will 
then decay into its renormalized ground state, giving up 



its excess binding energy U to the bulk electrons in the 
metal. The transition rate per unit time to deposit an 
electronic energy W to the metal is given by 

T(W) = (V'o, V5(E + W- H')VtJj ) . (52) 



From Eqs.(50) and (52) it follows that 



T(W)dW 



2nJJ z-(E + W) 



(53) 



Thus the binding energy U = — B is rigorously deter- 
mined by Eqs.(48) and (53) to be 



B 



T{W)dW 
2n J J W + B 

while Eqs.(51) and (53) imply 
P 



1 + zu 



where 



2tt 



In the perturbative limit, 



T{W)dW 
{W + B) 2 



lim (P/P c ) = lim (U/U c ) = 1 

h — *-oo h — *oo 



(54) 



(55) 



(56) 



(57) 



where the Casimir values Uc and Pc have been de- 
fined, respectively, in Eqs.(35) and (37). For the non- 
perturbative limit of atomic adsorption for a finite height 
h = h s , the exact results of Eqs.(54), (55) and (56) can 
be employed for making realistic estimates of the effective 
charge Z e ff\e\ of the atom. 



VI. NET CHARGE ON AN ADSORBED ATOM 

For estimating the effective charge on the adsorbed 
atom, we note the following: (i) Eq.(54) is an im- 
plicit equation for the binding energy B = f{B). 
(ii)_ If f = lim w ^ T(W) > 0, then f(B -► 0) « 
(%T i 1 2ir)ln(const / 'B) , where a constant "cut-off" must be 
placed in the logarithm, (iii) We employ the following 
simple dispersion formula for T(14 7 ): Over a bandwidth 
< W < E a , we consider T(W) = f to be uniform. 
Outside this interval, i.e. for electronic energies above 
the atomic energy cut-off E a = hto a in Table I, we con- 
sider r(W) to be negligible. Under this assumption, the 
binding energy Eq.(54) reads 



B 



— In 
2tt 



~B~ 



(58) 



if B « huj a . Similarly, Eq.(56) reads 



5 



\2ttB J 



(59) 



The probability of excitation P = Z e f f determines the 
effective charge via Eqs.(55), (58) and (59) 



Z eff ~ 



1 + ln(E a /B) 



(60) 



A plot of the effective charge versus the binding energy 
is shown in Fig. 2. 



0.25 



hand, the charge density functional approach is intuitive 
on a microscopic level. 

The evidence for Z e ff > is based on the change 
in metallic work function due to adding submonolayer 
adatoms, for thicker films the situation is complex. |17| ] 
The work function can increase as the film thickness in- 
creases for film layers above the monolayer. This experi- 
mental fact has (thus far) no simple theoretical explana- 
tion. 

APPENDIX 

In this appendix, the details of the mathematical 
derivation of results in Sec. V will be made explicit. One 
may define two projection operators Q and P obeying 



Q + P=l. 



(Al) 



0.2 - 



CD 0.15 



N 



The operator Q projects a wave function onto the un- 
perturbed ground state wave function ipo and may be 
written as 



I'M {H 



(A2) 



0.05 



0.005 



0.01 
(B/E 



0.015 



0.02 
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FIG. 2. The effective charge Z e ff\e\ on the atom is shown 
as a function of the binding energy B in units of the atomic 
excitation energy E a = hui a . See Table I. 



In terms of experimental binding energies and mean 
atomic excitation energies listed Table I, the effective 
charges can be computed numerically from Eq.(60). 



VII. CONCLUSIONS 

A dispersion theory has been presented for comput- 
ing the effective charge per noble gas atom adsorbed 
on metallic substrates. The general physical situation, 
shown in Fig.l, is complimentary to the charge den- 
sity functional approach for computing atomic binding. 
p8| , p9]| Both approaches yield a picture in which there is 
an atomic dipolc moment. The negative end of the dipolc 
moment is better viewed as the donation of a mean neg- 
ative electronic charge to the metal below. The positive 
end of the atomic dipole moment, is what has been called 
the net atomic charge Z e ff\e\. 

The advantage to the dispersion method is that it fits 
smoothly into the power law limit of the well known 
Casimir force. In the density functional approach, the 
"image charges" are added in "by hand". On the other 



If we define for any operator A, the projected operators 
A QQ = QAQ, A QP = QAP, A PQ = PAQ and A PP = 
PAP, then the operator may be written in a partitioned 
matrix form 



.4 



Aqq Aq P 
Apq App 



In particular, the resolvent operator 

I 



K{z) 



z-n 



may be written as 



V(A - ( n QQ( z ) n Qp( z ) 
K[Z) ~ {tZpq(z) Kpp(z) 



where 



G{z)l QQ = TZqq(z) = 1 QQ 



(z - £ n ) 



(A3) 



(M) 



(A5) 



(A6) 



From Eq.(A6), it follows that G(z) has a ground state 
pole at zq = £q with a residue given by \(^S> ni ipo)\ 2 ; i.e. 



G (z)-J?°'*»>j a „ : 



z . 



(z - Z ) 

Eqs.(A4) and (A5) imply 



z-T-Lqq -T~Lqp \{Kqq{z) Kqp(z) 
-TLpq z -Hpp ) \U P q(z) Upp(z) 



(A7) 








Ipp 



(A8) 



G 



from which 

(z - H QQ )TZ QQ (z) - H QP K PQ {z) = 1 QQ , (A9) 

and 

-HpqKqq(z) + {z- Hpp)K pq (z) = 0. (A10) 
Eqs.(A9) and (A10) imply 

(z - 7-Lqq - Hqp z _ npp n PQ^ n QQ( z ) = l QQ- 

{All) 

Finally, if H = H + V, (ip , V^o) = and W = H PP , 
then Eqs.(A6) and (All) read 

{z-E -(i, ,v{j-^}v^G(z) = l-, (A12) 

i.e. 

with a self energy 

(v^j^^J^o), 0414) 
which completes our derivations. 
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